Observation of sub-Bragg diffraction of waves in crystals 
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We investigate the diffraction conditions and associated formation of stopgaps for waves in crystals 
with different Bravais lattices. We identify a prominent stopgap in high-symmetry directions that 
occurs at a frequency below the ubiquitous first-order Bragg condition. This sub-Bragg diffraction 
condition is demonstrated by reflectance spectroscopy on two-dimensional photonic crystals with 
a centred rectangular lattice, revealing prominent diffraction peaks for both the sub-Bragg and 
first-order Bragg condition. These results have implications for wave propagation in 2 of the 5 
two-dimensional Bravais lattices and 7 out of 14 three-dimensional Bravais lattices, such as centred 
rectangular, triangular, hexagonal and body-centred cubic. 



The propagation and scattering of waves such as light, 
phonons and electrons are strongly affected by the peri- 
odicity of the surrounding structure [TJ [2J. Frequency 
gaps called stopgaps, emerge for which waves cannot 
propagate inside crystals due to Bragg diffraction. Bragg 
diffraction is important for crystallography using X-ray 
diffraction [3] and neutron scattering !4j. Diffraction de- 
termines electronic conduction of semiconductors [TJ [5J 
and of graphene [5], and broad gaps are fundamental for 
acoustic properties of phononic crystals [S] [7] and optical 
properties of photonic metamaterials [9] [10] . 

Bragg diffraction is described in reciprocal space by 
the Von Laue condition k ou t — k - ln = y, where fc ut, 
k j n are the outgoing and incident wave vectors and g 
is a reciprocal lattice vector. As a result, a plane exists 
in reciprocal space for which the Von Laue condition is 
satisfied, called a Bragg plane. When the incident and 
outgoing wave vectors are located on a Bragg plane these 
waves are hybridized, thereby opening up a stopgap at 
the Bragg condition. The boundary of the Brillouin zone 
is formed by intersecting Bragg planes and therefore gaps 
open on this boundary [ilj. When diffraction involves a 
single Bragg plane, we are dealing with well-known sim- 
ple Bragg diffraction, which corresponds in real space 
with the well-known Bragg condition: m\ = 2dcos(6). 
Here m is an integer, A is the wavelength inside the crys- 
tal, 6 is the angle of incidence with the normal to the 
lattice planes, and d is the spacing between the lattice 
planes. A stopgap is also formed when Bragg diffraction 
occurs on multiple Bragg planes simultaneously, which is 
called multiple-Bragg diffraction [TT], and is fundamental 
for bandgap formation [2j [12j [13] . Wave propagation in 
crystals is described along high-symmetry directions pQ. 
Multiple-Bragg diffraction has been recognized in high- 
symmetry directions at frequencies above the first-order 
simple Brag^g diffraction condition: m = 1, A = 2d or 
k out — — k in = g y ■ To our knowledge, multiple-Bragg 



diffraction has not yet been observed at frequencies below 
simple Bragg diffraction [T4"] . 

In this Letter we show that for high-symmetry direc- 
tions multiple-Bragg diffraction can occur at frequen- 
cies below the first order simple Bragg condition. As 
a demonstration we have investigated diffraction condi- 
tions for two-dimensional (2D) photonic crystals using 
reflectance spectroscopy. A broad stopgap is observed 
below the simple Bragg condition, depending on the sym- 
metry of the lattice. Our findings are not limited to light 
propagation, but apply for wave propagation in general, 
and therefore we anticipate similar diffraction for elec- 
trons in graphene [5], and sound in phononic crystals 

We have studied light propagation in 2D silicon pho- 
tonic crystals [T7]. Figure [IJa) shows a scanning electron 
microscope (SEM) image of one of these crystals from 
the top view. The centred rectangular unit cell has a 
long side a = 693 ±10 nm and a short side c = 488 ±11 
nm. The pores have a radius of r = 155 ± 10 nm and 
are approximately 6 fim deep. The photonic crystals are 
cleaved parallel to either the a-side or c-side of the unit 
cell. The cleavages define two directions of high sym- 
metry, TM' and TK, in the Brillouin zone, see Figure 
[IJ&). If light travels parallel with these directions, one 
expects simple Bragg diffraction from the lattice planes 
in real space (dashed lines in Figure [TJa)). A stopgap 
should appear that is seen in reflectivity as a diffraction 
peak. Because both directions are of high-symmetry, one 
naively expects that simple Bragg diffraction to give the 
lowest-frequency diffraction peak. 

We have identified the diffraction conditions of our 2D 
photonic crystals along the TM'- and Tif-directions us- 
ing reflectance spectroscopy [TS]. The photonic crystals 
arc illuminated with a supcrcontinuum white light source 
(Fianium SC-450-2). TE-polarized light is focused on the 
crystal using a gold-coated reflecting objective (Ealing 
X74) with a numerical aperture of 0.65, resulting in a 
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Figure 1: (color online) (o) SEM image of a 2D photonic 
crystal with a centred rectangular lattice. The white rectan- 
gle marks an unit cell with a = 693 ± 10 nm, c = 488 ±11 
nm and r = 155±10 nm. The blue and pink arrows mark 
two directions of high symmetry YK and FA/'. The blue and 
pink dashed lines mark real space lattice planes whose lowest- 
frequency simple Bragg diffraction occurs along the YK- and 
FM'-directions. (b) Reciprocal space of the centred rectangu- 
lar lattice (circles). The filled area is the first Brillouin zone, 
bi and 62 are primitive vectors. F, K, K' , M, and M' are 
points of high symmetry. The dashed lines are Bragg planes. 



spectrum angle-averaged over 0.447T ± 10% sr solid angle 
in air. By assuming an average refractive index (n=2.6), 
the angular spread inside the crystal is only 14°, corre- 
sponding to 0.067r± 10% sr solid angle. The diameter of 
the focused beam is estimated to be 2u>o = 1 A*m. Re- 
flected light is collected by the same objective, and the 
polarization is analyzed. The spectrum is resolved using 
Fourier transform infrared spectroscopy (BioRad FTS- 
6000) with an external InAs photodiode. The spectral 
resolution was 15 cm -1 , corresponding to about 10~ 3 rel- 
ative resolution. For calibration, spectra are normalized 
to the reflectance spectra of a gold mirror. 

Figure ^a) shows the bandstructure calculated us- 
ing a plane wave expansion method |19j and reflectiv- 
ity measured along the rM'-direction (black solid). The 
broad lowest-frequency measured reflectivity peak be- 
tween 4700 and 7300 cm -1 agrees well with the cal- 
culated stopgap. This reflectivity peak is caused by 
simple Bragg diffraction on the lattice planes indicated 
in the cartoon above, corresponding to the pink lattice 
planes in Figure [lja). One can also approximate the 
lowest-frequency simple Bragg diffraction condition from 
the dispersion with a constant effective refractive index 
(n e ff), obtained from the low- frequency limit |20j . This 
estimation is marked by the dashed vertical line and 
agrees well with the calculated stopgap. The two mea- 
sured peaks between 9800 and 11100 cm -1 agree well 
with a higher-frequency stopgap marked by a second blue 
area, caused by multiple-Bragg diffraction. The peaks 
appear at higher frequency than simple Bragg diffrac- 
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Figure 2: (color online) Measured (black, solid) and simulated 
(grey, dashed) reflectivity spectra, and calculated bandstruc- 
tures for TE-polarized light of a 2D photonic crystal along 
directions of high symmetry, (a) The measured and simu- 
lated lowest-frequency diffraction peaks in the FM'-direction 
match a calculated stopgap that occurs at the simple Bragg- 
diffraction condition. (6) The measured and simulated lowest- 
frequency diffraction peaks in the FTf -direction match a cal- 
culated stopgap and is caused by multiple-Bragg diffraction 
that occurs at a lower frequency than simple Bragg diffrac- 
tion. 



tion, as expected. The reflectivity of an incident plane 
wave on a finite size structure has been simulated with 
finite difference time domain (FDTD) simulations (grey, 
dashed). The agreement between the simulated and mea- 
sured reflectivity is gratifying. 

In Figure [2j&) we show the calculated bandstructure 
and measured reflectivity along the r indirection, where 
K is located on the edge of the Brillouin zone and B 
is located on the Bragg-plane between T and reciprocal 
lattice vector G\\. Two significant broad measured re- 
flectivity peaks are visible. The lowest-frequency peak 
between 5400 and 6900 cm -1 agrees well with a calcu- 
lated stopgap marked by the yellow area. This peak is 
caused by multiple-Bragg diffraction on the lattice planes 
indicated in the cartoon above the calculated stopgap, 
and is part of the two-dimensional band gap for TE- 
polarized light. The second reflectivity peak between 
8100 and 10000 cm -1 agrees with a second calculated 
stopgap (blue area). The flat bands in the dispersion re- 
lation, causing an impedance mismatch of coupling light 
into the crystal [TU1 [H] , likely broaden the observed peak 
(hatched area). This is supported by FDTD simulations 
of the reflectivity of an incident plane wave on a finite 
size structure (grey, dashed). The agreement between 
the simulated and measured reflectivity peak is gratify- 
ing. The measured peak is probably rounded-off as a 
result of the high-NA microscope objective. Note that 
bandstructure calculations and FDTD simulations ne- 
glect the dispersion of silicon. Scattering from surface 
imperfections becomes more important at higher frequen- 
cies, which could explain why the measured reflectivity 
peak is much lower near 10000 cm -1 . At any rate, the 
frequency ranges of the measured and simulated peaks 
agree very well. 

This second stopgap is caused by simple Bragg diffrac- 
tion on the lattice planes indicated in the cartoon above 
the calculated stopgap, corresponding to the blue lattice 
planes in Figure [TJa). The frequency of the simple Bragg 
diffraction condition based on an n e g is inaccurate be- 
cause a broad stopgap is already present at lower frequen- 
cies. The observation of a prominent diffraction peak 
caused by multiple-Bragg diffraction at a much lower fre- 
quency than simple Bragg diffraction is important. This 
result shows that even for high-symmetry directions such 
as the rif-direction, there can be a diffraction condition 
below simple Bragg diffraction, which we address as sub- 
Bragg diffraction [22] 

We have performed reflectivity measurements on pho- 
tonic crystals with a range of ~. Figure [3| shows the 
width of the diffraction peaks for the TM 1 - (a) and re- 
directions (b). The areas correspond with calculated 
stopgaps, such as in Figure [2j The dash-dotted line is 
the approximated frequency of lowest-frequency simple 
Bragg diffraction assuming a constant n e g. Note the very 
good agreement between the measured frequencies of the 
diffraction peaks and the calculated stopgaps. We ob- 
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Figure 3: (color online) Determined reduced width of the 
diffraction peaks (bars) and frequency of the maximum reflec- 
tivity (circles) for different ^. The filled areas are calculated 
stopgaps, color-coded as in Fig. [5] (a) Reduced frequency 
of the diffraction peaks for the rM'-direction. (6) Reduced 
frequency of the diffraction peaks for the TA'-direction. 



serve for the rif-direction diffraction always appearing 
at a lower frequency than simple Bragg diffraction. This 
observation confirms the robustness of sub-Bragg diffrac- 
tion. 

The existence of sub-Bragg diffraction can be explained 
by considering the lattice in reciprocal space, see Figure 
[TJfe). For the rX-direction we observe in reciprocal space 
two points of high symmetry: K and B. K is located on 
the Brillouin zone boundary, at the intersection of two 
Bragg planes corresponding to the von Laue conditions 
between T and Gio, T and G\\. Thus at K we have 
a multiple-Bragg diffraction condition on both Bragg 
planes. B is located at the Bragg plane (dashed line) 
that satisfies the von Laue condition between T = Goo 
and G\i resulting in simple Bragg diffraction. Since B 
is located outside the Brillouin zone, the simple Bragg 
condition occurs at higher frequency than the sub-Bragg 
condition. From this figure we describe three conditions 
for sub-Bragg diffraction: i. The diffraction condition 
corresponds to a point on a corner edge of the Brillouin 
zone, giving rise to multiple-Bragg diffraction, ii. The 
incident wavevector should be along a high symmetry di- 
rection, which is satisfied by considering only reciprocal 
lattice vectors Gkhi, for which \h\, \k\, \l\ < 1 or equiva- 
lent. Hi. Sub-Bragg diffraction can only occur at a lower 
frequency than the simple Bragg diffraction condition. 

Using these three conditions, it becomes evident that 
diffraction conditions for M and M' correspond to simple 
Bragg diffraction for Gio and G^j respectively. K' satis- 
fies criteria i and in, however, it does not satisfy criterion 
ii. This diffraction condition belongs to multiple-Bragg 
diffraction in a direction of lower symmetry, similar to the 
observation in Ref. [12] . Therefore, sub-Bragg diffraction 
is only observed at K. In this case we have measured 
the reflectivity of photonic crystals that strongly inter- 
act with light. For our crystals, we find that for - > 0.07 
a stopgap opens at K, and for - < 0.07 flat dispersion 
bands appear. 
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Figure 4: (color online) Normalized sub-Bragg diffraction 
condition (purple) and normalized simple Bragg diffraction 
condition (black dashed) as a function of -. The symbols 
mark the reflectivity peaks of Fig. [2j&), assuming identical 
n e ff for both diffraction conditions. Sub-Bragg diffraction is 
satisfied for - > 1. Labels f,2,3 refer to - shown in the 
insets (bottom). Inset I shows the reciprocal lattice (circles) 
for ft = 1, giving a square lattice and no sub-Bragg diffrac- 
tion. Inset 2 shows the reciprocal lattice (circles) for - — \f2, 
resulting in sub-Bragg diffraction at K. Inset 3 shows the 
reciprocal lattice (circles) for - = \/3, resulting in sub-Bragg 
diffraction at K and K' . 



Up to now we have considered a centred rectangu- 
lar lattice with long side a, short side c, and - = \/2. 
However, sub-Bragg diffraction can be expected for any 
f > 1 23 . To illustrate this we have made an analyti- 
cal model to explain the sub-Bragg diffraction frequency. 
We calculate |r_ftT| and \TB\ as a function of where the 
frequency of the sub-Bragg condition is proportional to 
^7 \TK | and the frequency of the simple Bragg condition 
is proportional to ^^jrS], where cq is the vacuum ve- 
locity. The results are shown in Figure |4j When 2 — » oo, 

sub-Bragg diffraction occurs at j^j = \. Inset 1 shows 
the reciprocal lattice for - = 1, corresponding to the 
square lattice. In this case \TK\ — \TB\ and therefore 
sub-Bragg diffraction and simple Bragg diffraction occur 
at the same frequency, violating condition Hi. Inset 2 
shows the reciprocal lattice for " = \/2, corresponding 
to the experimental conditions of the structures inves- 
tigated by us. Inset 3 shows the reciprocal lattice for 
- = \/3, corresponding to the triangular lattice. All 
three conditions for sub-Bragg diffraction at K are ful- 
filled. There is also a sub-Bragg diffraction condition for 
K' . It may seem that condition it is violated because the 
rif'-direction corresponds to G 2 y. However, because of 



the rotational symmetry of the Brillouin zone, K = K' 
and the diffraction conditions in the G 2 x-direction are 
identical to the Gn-direction, and therefore condition ii 
is satisfied. In a similar experiment performed by 24J a 
diffraction peak was observed at K. However, these ex- 
cellent experiments were compared with bandstructures 
between TK, since it was not recognized that there is 
also a diffraction condition at B. For the centred rectan- 
gular lattice, one must calculate bandstructures between 
TB to get accurately estimate the width of the stopgaps. 
This is evident from the bandstructures in Figure 
by comparing the width of the stopgaps when one would 
consider only TK instead of TB. 

In the case of three-dimensional (3D) crystals, if a Bra- 
vais lattice has a planar cross-section that can be de- 
scribed by a centred rectangular lattice along a direction 
of high symmetry, sub-Bragg diffraction will occur. For 
2D Bravais lattices sub-Bragg diffraction can occur for 
2 out of 5 Bravais lattices; centred rectangular and tri- 
angular (which is a special case of centred rectangular), 
see Figure |4| There are 7 out of 14 3D Bravais lattices 
that have a planar cross-section that can be described 
by a centred rectangular lattice in a direction of high 
symmetry; body-centred cubic, body-centred tetragonal, 
base-centred orthorhombic, body-centred orthorhombic, 
face-centred orthorhombic, base-centred monoclinic and 
hexagonal. We predict that sub-Bragg diffraction can 
occur for these 7 Bravais lattices. 

Sub-Bragg diffraction is valid for any kind of wave- 
propagation in structures that fulfill the symmetry con- 
ditions. Therefore we predict that for X-ray spectroscopy 
on crystals a sub-Bragg diffraction peak can be observed. 
As multiple-Bragg diffraction is required for photonic 
band gap formation, hence sub-Bragg diffraction can af- 
fect band gap formation [TS]- Indeed, the sub-Bragg 
diffraction condition is part of the 2D TE-band gap in tri- 
angular lattices [10]. For elastic wave diffraction a prop- 
agation gap is formed at the sub-Bragg condition and 
therefore also for phonons and for relativistic electrons, 
such as the case of graphene, which has a triangular lat- 
tice. 
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